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Abstract 

Let Ay = f, A is a linear operator in a Hilbert space H, y 1. N{A) := {u : 
Au = 0}, R{A) := {h : h = Au,u £ D{A)] is not closed, \\fs - f\\ < 5. Given fs, 
one wants to construct us such that lim^^g ll''^^ ~ y|| = 0- A version of the DSM 
(dynamical systems method) for finding us consists of solving the problem 

us{t) = -us{t) + T^J u(0) = no, (*) 

where T := A* A, Ta := T + al, and a = a{t) > 0, a{t) \ as t ^ oo is suitably 
chosen. It is proved that us := us{ts) has the property lim5_»o ||^<5 — uW = 0. Here 
the stopping time ts is defined by the discrepancy principle: 

e-^'-'^ais)\\Q-^\^fs\\ds = c6, (**) 

c G (1, 2) is a constant. Equation (*) defines ts uniquely and lim^^g ts = oo. 

Another version of the discrepancy principle is proved in this paper: let as be 
the solution to the equation a||Q~^/5|| = c6, and ts be the solution to the equation 
a{t) = as- If limt^oo = 0, then lim5_,o \\us — y\\ = 0, where us := us{ts) and 
us{t) solves (*). 



1 Introduction. 

Let A he a linear bounded operator in a Hilbert space H (or in a Banach space X), and 
equation 

Au = f (1) 
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be solvable, possibly non-uniquely. Let N{A) = N and R{A) denote the null-space and 
the range of A, respectively. Denote by y the (unique) minimal-norm solution to (Q), 
y ± N. Given fs, \\fs — f\\ < S, one wants to find a stable approximation us to y: 

lira \\us — y\\ = 0. (2) 

There are many ways to do this: variational regularization, quasisolutions, iterative 
regularization (see, e.g., PP, 0, |H])- 

Here we study a version of the dynamical systems method (DSM) for finding us: 

Mt) = -us{t) + T;il)A*fs, us{0) = uo, (3) 
where T := A* A is selfadjoint, Ta := T + al, I is the identity operator, 

< a{t)] a{t) \ as t ^ oo; lim - = 0, d := (4) 

t^oo a at 

The element us in ^ is us{ts), where us(t) is the solution to Q, and ts, the stopping 
time, is found from the following equation for the unknown t: 

f e-^'-'^a{s)\\Q-l^^fs\\ds = c6, c e (1,2), (5) 

^ 

where Q := AA* is selfadjoint, Qa = Q + al, c is a constant, and \\fs\\ > c5. This 
equation we call a discrepancy principle. About other versions of discrepancy principles 
see 0, ^-m. 

The main result of this paper is the following theorem. 
Theorem 1. Assume that (@)) holds and 

lime*a(t)||Q„Y,)MI = oo. (6) 

Then equation (0) has a unique solution ts, 

\\mts = oo, (7) 

and Q holds with us := us{ts)- 

Remark 1. Assumption (P) is always satisfied if fs ^ R{A). Indeed, 

a'^d{Exfs, fs) 



a— >0 a— >0 Jq 



\pfEr>o, 



where E\ is the resolution of the identity, corresponding to the selfadjoint operator Q, P 
is the orthoprojector onto the null spac e N{Q ) ofQ, N{Q) = N{AA*) = N{A*) := N* , 
and \\PN*f5\\ >0 if fs ^ R{A), because R{A) = (N*)^ . 



In Section 121 we prove Theorem^ and Theorem |2l which says that holds without 

a2{t) 



assumption (0) but with an extra assumption limj^oo -^jj) = 0. 



2 Proofs. 

Let h{t) := a{t)\\A-^^^fs\\ := a{t)g{t). 
Lemma 1. Assume ©■ Then 

e'hif) 



lim 

/o e'h{s)ds 

provided that 

lim ^ = 0. 

Proof of LemmaUi Apply L'Hospital's rule to (jH)) and conclude that (jHJ holds if 

e*/i(t) + e*/i 
lim -— = 1. 

t-»oo e*/i(t) 

Relation ^ holds if © holds. 
Lemma n is proved. 

Lemma 2. Relation © /ioMs ^/ ano? only if 



Proof of Lemma\^ We have 

, (h^y , /i 

lim = 2 lim — . 

t^oo /l^ t^oo h 

Lemma 121 is proved. 

We have h = ag. Therefore 

jh'y _ {a')'g' + a\g'r _ ja')' ^ jg')' 

a^g'^ g"^ 

If holds, then Lemma 2 implies 
Let us prove that 

t-+oo g-^ 



Let Eg be the resolution of the identity of the selfadjoint operator Q. Then 
and 

dp{s) \d\ a{t)dp{s) 



(gr = -2a 



[s + a{t)]^ a Jq [s + a{t)] 



\m f' dpjs) _ 2\d{t)\ , 
- a{t) J, [s + a{t)]^ a{t) ^ ^ ^' 



(17) 



where we have used the monotone decay of a{t), which imphes — d = |d|. Our results 
remain valid for b = oo and their proofs are essentially the same. 
From dH), (Uni) and dm) it follows that 

(a'^Y Idl 
< lim < 2 lim ^ = 0. (18) 

We have proved the following lemma. 
Lemma 3. // (@]) holds, then 

li,„<^ = (19) 

From Lemmas 1-3 we obtain the following result. 
Lemma 4. // (jj)) and (jH)) hold, then 

r'e-^'~''^h(s)ds , , 

lim ^-^ / ^ = 1. 20 

t^oo h{t) ^ ' 

Proof of Lemma^ Relation ()20|) is equivalent to 

. 1. (21, 



/o e'h{s)ds 
By Lemmas 1-3, relation (j?H) holds if 

lim - = 0. (22) 

This relation holds by assumption 

Lemma m is proved. □ 

Now let us prove Theorem ^ 
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^(a) := a^Q-'fsf = t = cl5\ (24) 



Proof of Theorem^ By Lemma E] equation for large t can be written as 

amQ:lt)fs\\ = + 0(1)], t ^ oo. (23) 

Denote c[l + o(l)] by Ci. If c G (1,2), then ci G (1,2) for sufficiently large t. Fix a 
sufficiently large t and denote a(t) := a. Consider the equation 

{a + s)2 

The function 'il){a) is defined on (0,oo), is monotonically growing and continuous, 

^(+oo) = [ dp= Wfsf > cl5\ ^-(+0) = \\Pf8\\\ 
Jo 

where P, P := Eq — Eq_, is the orthogonal projection onto the null-space N{Q) of the 
operator Q, N{Q) = N{AA*) = N{A*) := A^*. Since \\fs - f\\ < 6, and / G i?(A) ± 
N{A*), we have 

\\Pfs\\<\\Pifs-f)\\ + \\Pf\\<S, (25) 

because ||P|| < 1 and ||-P/|| = 0. Therefore equation has a unique solution a = as, 
lim^^o as = 0. Consequently, equation (0) has a unique solution ts, which can be found 
from the equation 

a{t) = as. (26) 

Since a(t) \ 0, we have 

limt^ = oo. (27) 

Finally, let us prove that the element us := us{ts) satisfies relation Q. First we give a 
simple proof of this relation under an additional assumption ()39|). see below. Then we 
give a more complicated proof which does not use any additional assumptions. 
We have 



Thus 



us{t) = e-'uo + I e~^'-''>T;^l^A*fsds. 

J 

\usit) - y\\ < e-'\\uo\\ + \\ j\~^'-^^T;^^^A* {fs - f)ds\\ 

+ II f e-^'-^^T-^l^A*fds - y\\ := j, + n + Js- 

'J 



(2J 



(29) 



Clearly, limj^oo Ji = 0. We have ||T„-^A*|| <^,a>Q. 

Indeed, A* = A*Q-^ = UQ^Q"^ , where U is an isometry, and we have used the 

formula T'^A* = A*Q~^ and the polar decomposition A* = UQ^. Thus ||T-M*|| < 

1 

WQ^Qa^W = sup5>o^ = 27s- Consequently, 

32 < — 30 
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Furthermore, / = Ay, so 

II j^e-^'-^^T-l^A*Ayds - y\\ < \\y\\e-' + l\-^'-^^a{s)\\T-l^y\\ds, (31) 



and 



limalr-^,f = lim / = = Q, (32) 

where Fg is the resolution of the identity corresponding to the selfadjoint operator T = 
A* A, V is the orthogonal projection onto N{T) = N{A), and Vy = because y -L N. 
It follows from (jSH) and ^ that 



t— >oo 

Indeed, if 



lim js = 0. (33) 



V9(s) := a{s)\\T^^l^y\\ as s oo, (34) 
then, as follows from Lemma 6 (see below), we have 



lim / e"(*"'V(s)rfs = 0. (35) 

If 

sr 

lim , = 0, (36) 

lim ts = oo, (37) 

then dSni), dSIIl) and (jSni) imply relation with us = us{ts). 
Equations and pHjl show that 

aWQa^fsW = ci6, l<ci<2, a:=as. (38) 

We do not have a proof of the relation (jHBj) . and it is not known, in general, if ()36p 
holds without additional assumptions. We can prove (IHUj) under the following additional 
assumption: 

y = T'^v, 0<7<^. (39) 
Lemma 5. // (|39|) holds and a = solves (|38p . t/ien (j36|) holds. 
Proof of Lemma{^ We have / = Ay = AT'^v, and (jHHj) implies 

Ci(5 



<V^(llQa-/(/-/^)ll + IIQa".'^3-"^ll) 



(40) 



where we have used the obvious estimate HQ^ ^|| < ^, which holds for any a > 0, and the 
commutation relation Q^^^ = ^T~^. Therefore 

(ci — 1)— — < max = c||^;||a7 — > as 5 — > 0, (41) 

^>o s + as 

where c = 0(7) > is a constant, Ci > 1, maXs>o — = c('~f)a'^~2^ and c = cH) = 
jf{l — py~'P with p := 7 + |. 

Lemma is proved. □ 

However, we want to prove relation (j2)) with us = us{ts) (where ts solves PU]) . and as 
solves (jHHj) ) without any extra assumptions of the type flH^ . 
To accomplish this, we propose the following argument. 
Equation (j38|) can be written as 

c,6 = \\aQ~'fs\\ = Ua + Q- Q)Q-'fs\\ = WQQ-'fs - fsW = \\AT-'A*fs - fs\\. (42) 

Let us denote ws := T~^^A* fs, where as is the solution to equation (|^. It is proved in 
p.22] that 

lirn \\ws — y\\ = 0. (43) 

\imus{ts)=y (44) 
0^0 



By formula (j28j) the relation 

holds due to the following lemma. 
Lemma 6. // 



lim q{t) = q, (45) 



then ^ 

lim / e-(*-")g(s)rfs = (46) 

Proof of Lemma\^ Take £ > arbitrarily small. Using ()45|). find t(£:) such that 

\q{t)-q\<^- for t > ^(e). 

Then 

/:=| /" e^^'~'^{s)dsq\<\ [ \'^'-'\{s)ds\ + \ [ e-^'''\{s)ds- q\ := h + h. (47) 

Jo Jo Jt{e) 

Take ^(e) sufficiently large, so that 

/i<| for t>r(£), (48) 
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and ^ ^ 

l2< [ e-^'-'^\q{s)-q\ds+ \ [ e'^'^Usq - q\ 

Jt{e) Jtie) (49) 

< |(1 - e-(*-*(^)) + |g|e-(*-*(^» < | for t > r(£). 

Then 

/<£ if t>r(£). (50) 

Lemma ini is proved. □ 

Theorem Q is proved. □ 

We can simphfy our arguments if the following extra assumption on a{t) is made in 
addition to (jD): 

lim ^ = 0. (51) 

Theorem 2. Assume (@)) and (j5ip . Let so/tJe (j38|) and ts solve (j26|) . T/ien relation 
((21) /loWs wt/i U5 := ^5(^5), where us{t) is given in 



Proof of Theorem\^ We have 

1^5(^5) - y\\ < \\u6its) - ws{ts)\\ + Wwsits) - y\\, (52) 
where ws{t) := T~^j-^A*fs. It is known (|^^, p. 22) that 

Let us prove that 
We have 



\im\\wsits)-y\\=0. (53) 
0— >o 



hm \\usits) -wits)\\ =0. 

(5—^0 



^^(t) = e"*Mo + / e~^*"')M;5(s)(is, (54) 
Jo 

and 

/ e-^'-'^W5is)ds = e-^'''^ws{s) - / e-^'-'^wsis)ds. (55) 

io "''0 

Prom (inD) and we get 

\\us{t)-nJs{t)\\<ce-'+ [ e-(*-^)||^i;,(s)||rfs. (56) 

Jo 

From (j5l|l and Lemma El it follows that 

lim \\us{ts) - ws{ts)\\ = 0, 

0— +0 



provided that 

lim = 0. (57) 

To prove we use assumption and the following estimate: 

where Fx is the resolution of the identity corresponding to the operator T (cf formula 
(32)). Theorem 121 is proved. □ 

Remark 2. Assumption ()5H) holds, for example, if a{t) = ^^^"^^^t for sufficiently large t, 
t > T, Co, Cib > are constants, b G (0, 1), and T > is an arbitrary large fixed number. 
On the initial finite interval t G [0,T] the function a{t) > may decay arbitrarily fast. 
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